r 

NASA Technical Memorandum 102140 

, Nf. f ; : • 

Frequency Domain Model for Analysis 
of Paralleled, Series-Output-Connected 
Mapham Inverters , 


Andrew S. Brush 
Sverdrup Technology, Inc. 

NASA Lewis Research Center Group 
Cleveland, Ohio 

Richard C. Sundberg 
General Dynamics Corporation 
San Diego, California 

and 

Robert M. Button r " 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

24th Intersociety Energy Conversion Engineering Conference 
cosponsored by the IEEE, AIAA, ANS, ASME, SAE, ACS, and AIChE 
Washington, D.C., August 6-11, 1989 


iw\s/\ 

{NASA-TM- 1021 40) FREQUENCY DOMAIN MODEL FOE 
ANALYSIS OF PARALLELED, 

SERIES-OUTPUT-CONNECTED MAPHAM INVERTEES 
/NASA. Lewis Research Center) 8 p CSCL 09A 

G3/33 


N89-261 49 

Unclas 

0220025 


I 



FREQUENCY DOMAIN MODEL FOR ANALYSIS 
OF PARALLELED, SERIES-OUTPUT-CONNECTED MAPHAM INVERTERS 


Andrew S. Brush 1 

Sverdrup Technology, Inc. 
NASA Lewis Research Center Group 
Cleveland, OH 44135 


Richard C. Sundberg 

General Dynamics Corporation 
P.O. Box 85990 
San Diego, CA 92138 


Robert M, Butt mi 

National Aeronautics m*! Space Administration 
Lewis Research Center 
Cleveland, OH 44135 


ABSTRACT 

The Mapham resonant inverter is characterized as a two- 
port network driven by a selected periodic voltage. The 
two-port model is then used to model a pair of Mapham 
inverters connected in series and employing phasor voltage 
regulation. It is shown that the model is useful for 
predicting power output in paralleled inverter units, and for 
predicting harmonic current output of inverter pairs, using 
standard power flow techniques. Some examples are 
compared to data obtained from testing hardware inverters. 

INTRODUCTION 

Power systems operating at a bus frequency of 20 kHz are 
currently under consideration for power distribution in 
spacecraft 1]. There are also potential applications in 
aircraft and other locations requiring compact 
electromechanical actuation systems[2j. Whatever the 
application, inverters will be required to convert DC source 
energy to 20 kHz. A leading inverter topology is the 
resonant circuit developed by Mapham[3], which offers high 
efficiency and low output distortion. In any such power 
system, multiple 20 kHz drivers will be employed for 
reliability and flexibility. Where multiple AC sources are 
connected in parallel to a system, the problem of regulating 
system voltage and power flow arises. A research program 
was undertaken at the NASA Lewis Research Center to 
characterize the performance of paralleled resonant 
inverters, to develop a model useful for control of inverters, 
and to increase our understanding of frequency domain 
analysis of resonant inverters. 

The functional block to be modeled was the Main Inverter 
Unit (MIU), containing two Mapham inverters connected in 
series and employing phasor voltage regulation. The results 
of modelling an MIU were compared to data obtained from 
the PMAD/PV testbed at LeRC[4], which contains three 
12 kW MIU’s. 


1 Work supported by Contract NAS3-25266, Task Order 8801, 
"PMAD/PV Testbed Support." 


MODEL DEVELOPMENT 

Time domain simulations such as EASY5[5] have been very 
successful in modelling the operation of resonant power 
converters[6][7J. However, a large model containing two 
MIU’s, a transmission line, and a load takes a large amount 
of computer CPU time to arrive at a steady state condition. 
In order to obtain a model which could be computed fast 
enough to be employed in on-line controllers, it was decided 
that a frequency-domain model would be employed. 

In standard references on AC utility power system analysis, 
such as Stevenson[8], a power source is typically modeled in 
the frequency domain as a voltage source in series with an 
impedance. This arrangement of source and impedance is 
similar to that described by Oruganti[9] as he proposed to 
model an inverter as a filtered voltage source. The 
advantages of taking the filter approach lie in the fact that 
it allows the use of standard load flow calculation 
techniques, and that the computations required to determine 
power transfer through the filter are simple. The issues to 
be resolved are selection of the proper filter parameters, 
and an appropriate input waveform. 

Filter 

The Mapham inverter power circuit to be modelled is shown 
in Fig. 1. This inverter is represented in the frequency 
domain by the equivalent circuit of Fig. 2. This is similar 
to models presented by Oruganti, and by Tsai and Lee[8], 
with the addition of the series capacitor and transformer 
reactances. Resistors have been inserted to model losses. 
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flowing in each mode. 
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Figure 2 Frequency Domain Circuit Model 


The main new step taken in this model is the selection of 
the value to use for the resonant inductance. Oruganti used 
L in his model, while Tsai and Lee used L/2 as an 
approximation. As observed by Sundberg, et. al.[ 10], neither 
of these values is accurate, since the effective resonant 
inductance of a Mapham inverter varies with the operating 
condition of the inverter. At no load, an accurate 
approximation of the effective resonant inductance L cff can 
be determined by measuring the period T^ of the inductor 
currents and noting that the resonant period is determined 
by the tank capacitor C r and the effective inductance: 


= 27tJl^c7 
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Unfortunately, when the inverter is loaded, the load 
becomes part of the resonant circuit, complicating the 
calculation of resonant period, and the above equations do 
not apply. 

To understand the mechanism which determines L cff , 
consider the inductor current waveforms of Fig. 3. During 
the time periods labelled Ml, current flows in only one path 
between the DC supply terminals, through legs A and D, or 
through legs B and C. The series inductance of path AD or 
BC is equal to L. During the periods labeled M2, all four 
inductors conduct, creating two parallel paths, both of 
inductance L. In this mode, the resonant inductance is 
clearly L/2. We propose to determine L cff as an average of 
L and L/2, weighted by the proportion of the total i 2 t 



Consider a half period starting at the beginning of the mode 
M2 in at time 0, transitioning into mode Ml in at time M, 
and continuing to time T/2. The effective inductance will 
be found by: 



+ i 

2 (^M2 + ^Ml) (^M2 + ^mJ 


I M2 is the total i 2 t flowing when all inductors are conducting, 
or when effective inductance is L/2. I M1 corresponds to an 
effective inductance of L. Since each half period is 
symmetrical, there is no need to average over a full period. 

Figure 4 shows a comparison between no load L eff 
calculated by the method of equations (1) and (2) with that 
found by (3), (4), and (5) at various switching frequencies, 
for the inverters in the testbed. It is seen that both methods 
arrive at a no load L eff of 0.86L when the switching 
frequency is 20kHz. 

In order to determine values of L C ff to use in the model, 
inductor current data from EASY5 simulations of a testbed 
inverter was used to produce the map of Lcff vs. load given 
in Fig. 5. Assuming that inverter output voltage is 235 V, 
data from the map is approximately fitted by the equation, 

<-eff * L, + I I out lysing (6) 

where 7 = 0.006 for the testbed inverter, and 0 is the angle 
by which inverter output voltage leads output current. 

With all circuit values determined by analysis or 
measurement, the circuit of Fig. 2 is represented by the two- 
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Figure 4 Effective Resonant Inductance vs. Normalized Switching 
Frequency 


Figure 3 Inductor Currents 
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Figure 5 Resonant Inductance vs. Load for various power factors. 


port of Fig. 6. Standard circuit analysis techniques are used 
to calculate the complex hybrid circuit parameters A, B, C, 
and D of the model two-port, such that: 

E in - AV out + BI out (7) 

Ii„ ■ CV out + 0I out (8) 

Of course, the parameters A, B,C, and D are dependent on 
frequency, and need to be calculated for each harmonic. 


Driving Voltage 

Unlike the two-port parameters, which represent the actual 
characteristics of the inverter circuit, the input voltage E in 
is a made-up quantity which cannot be measured during 
inverter operation. Previous authors have suggested 
selecting a driving voltage based upon the operation of the 
inverter. Oruganti used the square wave of Fig. 7(a), which 
is suitable for analysis at the fundamental, but is a rough 
approximation at other harmonics. The wave of Fig. 7(b) 
was calculated by Tsai and Lee from simulation data, based 
upon a model with fixed Leff of L/2, from the equation, 


d (i 


AD 


*BC) 
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+ v r 


(9) 


dt 


where v c is the voltage across the tank capacitor. Clearly, 
any change in the inverter operating condition will change 
the inductor currents, thereby changing the calculated e in . 
Also, as previously observed here, L cff changes with 
operating condition. 

It seems reasonable that a map of harmonic voltages and 
phase angles could be produced for all inverter loading 



Figure 6 Two-port representation of a Mapham inverter. 



conditions using (7) and the simulation data used to produce 
Fig. 5. Although this was rejected for the present model in 
favor of an E in of fixed harmonic composition, we are 
proceeding to evaluate the advantages of such a map. The 
E^ used in this model was calculated from actual measured 
open-circuit terminal voltage of an inve/ter ijj the testbed. 
For each harmonic considered, E in = AV^. V out is a vector 
determined from measurements of open circuit voltage 
made using a harmonic analyzing voltmeter. The e in (t) wave 
used in this model (Fig. 7(c)) is the sum of the voltages at 
the fundamental, and harmonics 2, 3, 5, 7, and 9. 

It has been observed that the output current of an inverter 
is linked to the phase relationship between the output 
voltage and the clock which times the firing of SCR’s in the 
inverter. In this model, we tested the theory that the phase 
of the fundamental of E in is analogous to the clock phase. 
While the phase of each harmonic in E in remains constant 
with respect to the fundamental, the fundamental is 
considered to be in phase with the inverter clock, which is 
variable in phase. This is not an exact assumption, but was 
deemed adequate for a simple model. 

Assembly of MIU model 

The testbed MIU’s consist of two complete inverters such 
as that in Fig. 1 connected with the transformer outputs in 
series. To model a MIU, two two-port networks 
representing the circuit of Fig. 2 are connected as shown in 
Fig. 8. The transformer connections impose the constraint 
that the output current of each inverter is identical to the 
MIU output current. When the two-port equations are 
combined with this constraint, MIU output current is given 
by: 


*miu 




AjBj + AjBj 


(10) 
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Figure 8 Model of Complete MIU. 


For given conditions of line voltage V linc and clock phases 
(which determine the phase of Ej and E 2 ), the 
corresponding I MIU must be found iteratively. This is due to 
the fact that the inverter operating conditions, and thence 
the ABCD parameters, are determined by the inverter 
output. A computer subroutine has been written which 
performs all the required calculations to find MIU output 
current at the fundamental, and at harmonics 2, 3, 5, 7, and 
9. This subroutine has been used in a loadflow of the 
simple system shown in Fig. 9. The loadflow used standard 
linear iterative methods to identify the line voltage and 
angle which would achieve conservation of real and reactive 
power. The following section compares results from the 
loadflow to data obtained from the testbed. 

RESULTS 

The loadflow of Fig. 9 was run using input data 
corresponding to actual testbed runs for which data was 
available. The usefulness of this model in predicting power 
output and harmonic current flow was checked against the 
hardware data. In the following discussion of results, the 
term "clock phase" will be used to refer to the phase angle 
of the clock of inverters #1 #3, and to the phase of the 
fundamental of E t and E 3 in the model. "Regulation angle" 
will refer to the angle by which the clocks of the inverters 
#2 and #4 lag inverters #1 and #3, respectively, and to the 
angle by which the fundamental of E 2 and E 4 lag E t and E 3 , 
respectively (Fig. 10). 
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Figure 9 System Model. 


MIU "B" 


Reference Angle 
Inv. *1 Clock (El) 
lnv. *3 Clock (E3) 

Regulation Angle A , nv . « 2 Clock (E2) 
Regulation Angle B 

Inv. *4 Clock (E4) 

Figure 10 Phasor Relationship between inverter clocks. 

Power Sharing * Phase Control 

In this test, the regulation angle of each MIU was 
maintained at zero, while the difference between the clock 
phases of MIU "A" and "B" was varied. The results are 
presented in Figure 11. These plots show that the model 
has qualitative agreement with hardware data. The model 
is more sensitive to changes in phase angle than the 
hardware. Also, the PMAD/PV MIU’s attained exact 
power division at a relative phase angle slightly off from the 
expected zero. This is no doubt due to differences in the 
actual component values between MIU's, and differences in 
propagation delays in the firing and logic circuits of the 
MIU’s. 




Difference Between Clock Phases (Degrees) 
Figure 11 MIU Output vs. Clock Phase Difference 
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Power Sharing • Regulation Angle Control 

This test was accomplished by leaving the clock phase of 
both MItPs at zero, and varying the regulation angle of each 
inverter. Figure 12 shows that, as in the phase control test, 
the PMAD/PV MIU’s achieve even power division at an 
offset from zero, probably due to the same causes. At lower 
power outputs, the model corresponds reasonably well with 
the breadboard data, with the exception that, again, the 
model is more sensitive to changes in angle. At higher 
power output, the model diverges from the test data. 

In the breadboard tests, and in the loadflow, differences in 
regulation angle in excess of fourteen degrees were achieved 
by setting MIU "A" to maximum voltage output (zero 
regulation angle), and increasing the regulation of MIU "B" 
to the desired point. Above this point, MIU "B" is supplying 
inductive VAR, and is exceeding its current rating. In this 
case, the assumption that E remains fixed in phase with 
respect to the inverter clock seems to break down badly. 


Harmonic Loadflow 

Figure 13 shows the output of the original loadflow model 
when operating at a phase angle of five degrees, and no 
regulation. Note that the fifth harmonic current is 
significant, although the line voltage remains undistorted. 
The oscilloscope photograph of Fig. 14, taken under similar 
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Difference in Regulation Angle (Degrees) 
Figure 12 MIU Output vs. Regulation Angle Difference. 
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Figure 13 Unfiltered MIU Output Calculated by Harmonic Loadflow 

conditions, reveals that the prediction of fifth harmonic 
circulating currents is accurate. 

Since it is preferable to eliminate circulating harmonic 
currents, some testing was conducted in the testbed to 



evaluate the efficacy of applying series-resonant LC filters 
to the inverter output. Fig. 15 shows the result given by the 
loadflow when a filter composed of a 2 iiF capacitor in 
series with a 31 mH inductor, with a 10 Ohm resistor 
shunting the LC combination, is applied to each MIU 
output. Clearly, the filter reduced the harmonic currents. 
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Figure 15 Filtered MIU Output Calculated by Harmonic Loadflow 
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Figure 16 Filtered M 1 U Current From Testbed. 

Fig. 16 gives the MIU outputs given by the testbed with the 
same filters. The model results are very similar to the 
actual measured data, despite the fact that losses in the 
filter were not modeled. 

Note that the operating condition used for this test, with 
one MIU supplying no power, is unrealistic, and was 
selected for demonstration purposes only. In more normal 
conditions, with inverters evenly dividing power output, 
circulating harmonics are minimal without filters. 

MODEL PERFORMANCE 

In the comparisons discussed above, the model gave 
acceptable qualitative results. The model accurately 
represented trends in inverter output with respect to the 
inverter clock signals, although overestimating the sensitivity 
of the system to changes in clock phase angles. A good 
forecast of circulating harmonic current was obtained, and 
the performance of a model filter is comparable to that of 
a real filter. 

The deficiencies in quantitative performance may stem from 
oversimplification in selection of the model input voltage E. 
It is obvious that the input voltage of equation (9) will vary 
in harmonic content as the resonant period changes with 
load. Not so evident is that the input voltage E seems to 
vary in phase with respect to the actual inverter clock as the 
resonant period varies with load. Further investigations 
using hardware data in conjunction with EASY5 time 
domain simulations are underway, and should increase 
understanding of how to model the voltage E. 


In a frequency domain model of a Mapham inverter, the 
phase of the model input voltage E in is approximately 
analogous to the phase of the inverter clock. 

The impedance model gave reasonable performance at 
harmonics of the fundamental switching frequency, and 
should be useful as a tool for filter design. 

The method of calculating the effective resonant inductance 
L cff gives better results in frequency domain analysis than 
either L or L/2 used previously. 
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CONCLUSIONS 

The frequency domain model of MIU’s operating in parallel 
verified the idea that impedance-based techniques for 
analysis of real and reactive power transfer can be applied 
to resonant inverters. 

The circuit analysis and loadflow techniques found in any 

introductory text on power system analysis are applicable to ORIGINAL PAGE 

high frequency systems employing resonant inverters. 8LACK AND WHITE PHOTOGRAPH 
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16. Abstract 

The Mapham resonant inverter is characterized as a two-port network driven by a selected periodic voltage. The 
two-port model is then used to model a pair of Mapham inverters connected in series and employing phasor 
voltage regulation. It is shown that the model is useful for predicting power output in paralleled inverter units, 
and for predicting harmonic current output of inverter pairs, using standard power flow techniques. Some 
examples are compared to data obtained from testing hardware inverters. 
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